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ABSTRACT 

The recent discovery of a correlation between nuclear black hole mass, M,, and the stellar velocity 
dispersion, (7*(Gebhardt et al. 2000a; Ferrarese & Merritt 2000), in elliptical galaxies and spiral bulges, 
has raised the question whether such a relationship exists for AGN. Estimates of M, for many AGN, 
made using reverberation mapping techniques, allow exploration of the relationship between black hole 
mass, the host galaxy and the energetics of nuclear emission. However, since only a few AGN have both 
M, and c* measurements, we use the [OIII] A 5007 emission line widths on the assumption that for most 
AGN the forbidden line kinematics are dominated by virial motion in the host galaxy bulge. We find 
that a relation does exist between M, and [OIII] line width for AGN which is similar to the one found 
by Gebhardt et al. (2000a), although with more scatter as expected if secondary influences on the gas 
kinematics are also present. Our conclusion is that both active and inactive galaxies follow the same 
relationship between black hole mass and bulge gravitational potential. We find no compelling evidence 
for systematic differences in the mass estimates from reverberation mapping and stellar dynamics. We 
also find that for radio quiet AGN the radio power and black hole mass are highly correlated linking 
emission on scales of kiloparsecs with the nuclear energy source. 

Subject headings: galaxies: nuclei — galaxies: Seyfert — quasars: general — galaxies: kinematics and 
dynamics 



1. INTRODUCTION 

In the standard model for an active galactic nucleus 
(AGN), radiant energy is released through accretion onto 
a supermassive (10^ — 10^ Mq) black hole. However, until 
recently, there have been few mass estimates for nuclear 
black holes in AGN. Kinematic studies of nuclear gas disks 
{e.g. Harms et al. 1994) and stellar dynamical studies, 
using high spatial resolution ground-based spectroscopy 
{e.g. Kormendy 1988) and more recently with the Space 
Telescope Imaging Spectrograph (STIS) {e.g. Bower et 
al. 2000), have produced numerous central mass estimates 
mostly for normal galaxies. However, in application to 
AGN these techniques are severely limited. Since roughly 
one galaxy in a hundred has a luminous active nucleus even 
the nearest examples are relatively distant reducing the 
effective spatial resolution. Furthermore, stellar dynami- 
cal techniques require high signal-to-noise measurements 
of stellar absorption features which are often lost in the 
glare of a bright active nucleus. Also emission lines from 
ionized gas at small radii can be influenced by relativistic 
jets or winds associated with the AGN. 

These difficulties can be overcome, however, using re- 
verberation mapping techniques {e.g. Netzer & Peterson 
1997). Monitoring the variability of the continuum and 
emission line fluxes in broad line AGN, shows that the two 
are highly correlated but that the emission lines lag behind 
the continuum by times ranging from days to months. The 
lag is best explained as a light travel time effect and can be 
used to estimate the size of the Broad Line Region (BLR). 
If one further assumes that the BLR kinematics are largely 
Keplerian, one can use the line widths as a characteristic 
velocity and in turn determine the total mass contained 
within the BLR. Strong evidence for Keplerian motion ex- 



ists for at least one well-studied Seyfert, NGC 5548. Re- 
sults for several different emission lines, each with different 
lags and therefore spanning a range of distances from the 
nucleus, give consistent values for the central mass (Pe- 
terson & Wandel 1999). Reverberation studies have now 
yielded M, estimates for numerous Seyfert 1 galaxies and 
quasars {e.g. Kaspi et al. 2000). However, there have 
been some claims that the masses determined from these 
variability studies are systematically low (Wandel et al. 
1999; Ho 1999). 

Recent work has focused on the distributions of black 
hole properties and the relationship to their host galaxies. 
Kormendy & Richstone (1995) first found a correlation 
between M, and the mass of the spheroidal component, 
Mbui . In a comprehensive study using ground-based spec- 
troscopy and HST imaging, Magorrian et al. (1998), pre- 
sented strong evidence in support of the the M, ~ M^ui re- 
lation. More recently, two studies (Gebhardt et al. 2000a; 
Ferrarese & Merritt 2000) reported relationships between 
M, and (Te, the stellar velocity dispersion obtained within 
a large aperture extending to the galaxy effective radius, 
thus insensitive to the influence of the black hole. The 
correlation between these two parameters is remarkably 
strong suggesting a link between the formation of the bulge 
and the development of the black hole. 

In this paper we combine reverberation mapping mea- 
surements of M, in AGN with Narrow Line Region (NLR) 
gas and bulge stellar kinematic measurements correspond- 
ing to motion well beyond influence of the black hole to 
compare the nuclei of active and normal galaxies. We also 
investigate the possibility of a correlation between black 
hole mass and the radio luminosity in AGN. Section 2 de- 
scribes the data collected from the literature. Section 3 
presents the M,- a^ relation for AGN. In section 4 we 



consider the relationship between M, and the radio lumi- 
nosity and in Section 5 we summarize our conclusions. 

2. DATA 

We start with all AGN with black hole masses, mea- 
sured by the reverberation mapping technique: 17 Seyfert 
1 galaxies (as compiled by Wandel et al. 1999) and 17 PG 
quasars (Kaspi et al. 2000). For most of the galaxies, we 
use the black hole masses tabulated by Kaspi et al. (2000) 
evaluated using the mean H/3 profile width. M, values for 
three galaxies (Mkn 279, NGC 3516, and NGC 4593) were 
taken from Ho (1999). For as many objects as possible 
we obtained published measurements of cr*, FWHM[oiii] 
, and the radio flux. The adopted values are presented 
in Table 3. For the Seyferts, FWHM[oiii] is taken from 
Whittle (1992a). For the quasars FWHMjoiii] values were 
only adopted if the spectroscopy was of medium to high 
resolution (i? > 1500). cr* in AGN host galaxies has been 
determined using spectral regions around both the Ca II 
triplet and Mg6. The values used here have been published 
in Smith, Heckman & lUingworth (1990), Terlevich, Diaz 
& Terlevich (1990) and Nelson & Whittle (1995) (see Nel- 
son & Whittle 1995, for references on individual galaxies). 
Radio luminosities arc from Whittle (1992a) for Seyferts 
and from KcUerman et al. (1989) for PG quasars. We have 
scaled observations at 5Ghz assuming a power-law index of 
0.7 and assume Hq = 80kms~^Mpc~-'^. We also calculate 
i?, the ratio of the radio to optical flux, using the mean 
continuum luminosities from Kaspi et al. (2000). 

We point out that the stellar and gas kinematic measure- 
ments are from small aperture observations typically cov- 
ering 2 — 3" somewhat smaller than the apertures used to 
measure (Tg. Nevertheless, we can be confident that these 
measurements are unaffected by the presence of the black 
hole. First, since the AGN are at larger distances than 
the normal galaxies for which M, has been determined, 
the aperture typically corresponds to physical scales of a 
kpc or more, actually comparable to the bulge effective ra- 
dius. Second Nelson & Whittle (1996) found that Seyferts 
were offset from the Faber- Jackson relation {M^ui oc logcr*) 
for normal galaxies in the sense of having lower velocity 
dispersions than normal galaxies of the same Mf,„; . Their 
interpretation was that Seyfert bulges have lower M /L ra- 
tios than normal spirals, exactly the opposite of what one 
would expect if the stellar kinematics were strongly influ- 
enced by a massive nuclear black hole. 

3. THE M,- CT* RELATION FOR AGN 

Recently, Gebhardt et al. (2000b) included Seyferts in 
their plot of a^^ vs. M, (in some cases choosing different 
values for a^, than in Table 3) and conclude that AGN show 
no significant difference from the overall relation. Unfor- 
tunately, cr* values are available in the literature for only 
a hand-full of the AGN with reverberation mapping esti- 
mates for M, . Therefore a definitive comparison of AGN 
and normal galaxies is not possible with these data. How- 
ever, Nelson & Whittle (1996) have shown that for the 
majority of Seyfert galaxies, a moderately strong correla- 
tion between FWHM[oiii] and ct* exists indicating roughly 
equal absorption and emission line widths. Thus the [OIII] 
A5007 profiles are dominated by virial motion in the bulge 
potential. As might be expected, a fair amount of real (i.e. 



not due to measurement error) scatter does exist. Nelson 
& Whittle (1996) used the deviation from purely virial gas 
motion to investigate possible secondary influences on the 
NLR kinematics. Their results conflrmed conclusions from 
previous studies {e.g. Whittle 1992b, c) that the interac- 
tion of NLR gas with a relatively strong kiloparsec-scale 
linear radio source can produce non- virial gas acceleration. 
They also found a weak tendency for interacting systems 
and mergers to have broader emission lines than expected 
from purely virial motion, an issue that may be important 
in quasars which are more likely to be interacting systems 
(McLeod & Rieke 1995). Thus, keeping the issue of sec- 
ondary influences on the NLR kinematics in mind, we can 
proceed with the idea that the primary influence on the 
forbidden line widths in AGN is the bulge potential. 

Tabic 1 
AGN with Measured Blaek Hole Masses 
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[0111] FWHM References: 1 Appenzeller & Wagner (1991), 2 Brotherton 
et al. (1996), 3 Gelderman & Whittle (1996), 4 Heckman, Miley & Green 
(1984), 5 Marziani et al. (1996), 6 NW95, 7 Vrtilek & Carleton (1985) 
(1985), 8 Whittle (1992a) 



In Figure 1 we plot FWHMfom] /2-35 or cr, vs. M,. 
Filled symbols show the FWHM[oiii] values for the galax- 
ies in Table 3, • for Seyferts, ■ for the PG quasars, -|- 
signs show AGN for which a^, measurements exist and 
triangles show M, and Ue from Gebhardt et al. (2000a). 
The similarity in the trends between the data from Ta- 
ble 1 and the Gebhardt et al. (2000a) sample is striking. 
The long solid line is the relation determined by Gebhardt 
et al. (2000a) and the shorter dashed line is our fit to 



the AGN, using the ordinary least-squares bisector (Isobe 
et al. 1990) which gives good results for large uncertain- 
ties in both variables. The correlation for the AGN is 
moderately strong {R = 0.51, P{null) = 0.5%) with a 
slope 3.7 ± 0.7 and intercept -0.5 ± 0.1. The AGN re- 
lation is slightly shifted to larger widths for a given M, 
than the Gebhardt et al. (2000a) relation and is more 
scattered. Also, there seems to be no statistical differ- 
ence between the distributions of Seyferts and quasars. 
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Fig. 1. — The data from Gebhardt et al. (2000a) is shown as A 
and the soUd Hnc is the M, - Us relation derived there. Objects 
with (Tt from Nelson & Whittle (1996) are shown as +. Objects for 
which the [OIII] line widths have been converted to a are shown as 
• for Seyferts and as ■ and the dashed line is the fit. 

The slope on the fit agrees well with that found by Geb- 
hardt et al. (2000a) who derived a different value than Fer- 
rarese & Merritt (2000) (3.75 as opposed to 4.8). Geb- 
hardt et al. (2000a) discuss various explanations for the 
difference in the two values. Although our result supports 
a shallower slope, an analysis using stellar velocity disper- 
sion instead of emission line width would be more conclu- 
sive. 

We can be more quantitative about the similarity of the 
relation for elliptical galaxies and that for the AGN, by 
considering the origin of the scatter. The mean 1 — a error 
bar in log M,, plotted in the upper left, is -1-0.17/ — 0.20 
while the r.m.s. deviation of the filled symbols in Figure 
1, relative to the Gebhardt et al. (2000a) relation is 0.55. 
Since we expect that emission line widths can be deter- 
mined to accuracies of ~ 10% for the spectral resolution of 
the observations, we expect that a large component of the 
scatter is not due to measurement error but reflects real de- 
viations perhaps indicating secondary influences on NLR 
kinematics. If we calculate the difference from the Geb- 
hardt et al. (2000a) relation, this time in log FWHMjoiii] , 
we find a small offset, Alog FWHM[oiii] — 0.06, with a dis- 
tribution of r.m.s. width 0.14. Nelson & Whittle (1996) 
used differences in emission and absorption line widths, 
AW = logFWHM[oiii] — log2.35cr*, to investigate non- 
virial infiuences on emission line kinematics. For their en- 
tire Seyfert sample they found AW ~ 0.0 and an r.m.s. 
width of 0.2. Excluding interacting galaxies and objects 
with relatively luminous linear radio sources they found 



AW = —0.1, thus slightly narrower emission lines than 
absorption lines, and an r.m.s. width of 0.13. Thus it 
is reasonable that the scatter and the small offset indicate 
non-virial contributions to the line widths in some objects. 
The fact that the r.m.s. deviation is comparable to that 
found in Nelson & Whittle (1996) is consistent with this in- 
terpretation. Unfortunately, the current sample includes 
no Seyferts with strong linear radio sources and, due to 
their larger distances, the radio maps for the quasars are 
mostly unresolved on kpc scales (see e.g. Kellerman et al. 
1994; Kukula 1998). Therefore, no strong statements re- 
garding non-virial gas kinematics can be made. 

The relationship between M, and the FWHM[oiii] in 
AGN is precisely what we would expect if AGN follow the 
M, - (7, relation and if the primary influence on NLR 
kinematics is the bulge gravitational potential. Although 
one might point out that the AGN plotted with ct* do 
seem to lie systematically below the relation, the reality 
of such a shift based on only 7 objects is not compelling. 
Our primary conclusion is that the relationship between 
the mass of the central dark object and the depth of the 
gravitational potential in AGN is the same as for normal 
galaxies. Clearly more velocity dispersion measurements 
for AGN are needed. However, in fainter objects detailed 
studies of the [OIII] emission line proflles accounting for 
non-virial kinematic components could provide a valuable 
substitute. 
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Fig. 2. — Log Li4,i5MHz vs. M, for AGN. 

These results seem to conflict with previous reports that 
AGN have lower black hole to bulge mass ratios than nor- 
mal galaxies (Wandel et al. 1999) or alternatively that the 
reverberation mapping estimates of M, are systematically 
low (Ho 1999). Figure 1 suggests that neither of these 
is correct and that the differences must lie in the values 
of Mhui- A number of factors may have contributed to 
this apparent disagreement. First, the M, - (jg relation 
derived by Gebhardt et al. (2000a) has been evaluated us- 
ing only stellar dynamical black hole mass estimates from 
three-integral modeling. The result is that the black hole 
masses are about a factor of 3 lower than the ones pub- 
lished by Magorrian et al. (1998) which were used in Wan- 
del et al. (1999). Second Nelson & Whittle (1996) found 
that Seyferts are offset from the Faber- Jackson relation 



for normal galaxies having brighter bulges by on average 
about 0™6 at the same velocity dispersion. More gener- 
ally bulge-disk decomposition is a difficult task even for 
normal, nearby galaxies. The problem is compounded in 
Seyfert 1 galaxies by their smaller apparent bulge sizes, 
due to larger mean distances, and the bright point source 
component. Thus, much of the discrepancy can be resolved 
by acknowledging that as much or greater uncertainty ex- 
ists in the bulge luminosity of Seyferts as exists in the 
estimates of M, by any of the available techniques. 

4. RADIO LUMINOSITY AND BLACK HOLE MASS IN AGN 

Using reliable measurements of black hole mass in AGN 
we can compare parameters related to the energetics of 
AGN with the size of the central engine. For example 
Kaspi et al. (2000) found a strong correlation between the 
black hole mass and the mean optical continuum lumi- 
nosity, with the dependence M oc L°-^, suggesting that 
more luminous AGN radiate at a larger fraction of the 
Eddington luminosity. Nelson & Whittle (1996) found a 
correlation between ct* and the radio luminosity at 1415 
MHz, Li4i5, in Seyfert galaxies. This relation was linked 
to previously known correlations between radio power and 
absolute blue luminosity (Meurs & Wilson 1984; Edel- 
son 1987; Whittle 1992c) and between radio power and 
[OIII] hnc width (e.g. Whittle 1985). They speculated 
that the origin of such a correlation might result from 
the correlation of M, with bulge mass and that the ra- 
dio sources associated with larger central engines might 
be correspondingly more luminous. In Figure 2 we plot 
LiAibMHz vs. M,, plotting different symbols for radio- 
loud {R > 8.0,+), radio-intermediate(1.0 < R < 8.0, o) 
and radio-quiet AGN(i? < 1.0, •). The criteria are cho- 



sen rather arbitrarily, but nevertheless demonstrate some 
interesting differences. The solid line shows a fit to the 
most radio-quiet group. The correlation is quite strong 
{R = 0.82, P(null)=0.003%) with a slope 0.84 ± 0.12. A 
fit including radio-intermediate data (not shown) is some- 
what weaker (R = 0.58, P(null)=0.1%) and a fit to the 
entire sample weaker still (i? = 0.40, P(null)=2%). The 
results suggest that at least in radio-quiet AGN the radio 
luminosity, emitted on kpc scales is linked to the mass of 
the nuclear black hole. Nelson & Whittle (1996) noted 
other explanations for their result, including higher am- 
bient pressure in larger bulges leading to increased radio 
emissivity. We cannot distinguish between these possibil- 
ities and it is possible that both effects work together to 
produce the observed correlation. 

5. SUMMARY 

Using the assumption that NLR kinematics are predom- 
inantly due to virial motion in the host galaxy potential, 
wc find that AGN follow a relation between M, and (Te 
similar to normal galaxies. Although the AGN black hole 
masses are determined using reverberation mapping, we 
find no significant differences with results from studies us- 
ing stellar dynamical techniques. Satisfied that the black 
hole masses in AGN are reliable we can begin to examine 
the relationship of black hole mass with other properties 
of AGN. As an example, we have found that the radio lu- 
minosity in radio-quiet AGN is correlated with black hole 
mass relating parameters determined on vastly different 
size scales. 

I would like to thank Donna Weistrop for her support 
and comments on a draft of this letter. This work was 
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